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Abstract 


This report contains the results of research activities conducted under 
Grant NAG 5-389 during 15 May 1986 to 14 February 1987. The first 
chapter deals with spatial and temporal variations of surface radiative 
temperatures of the burned and unbumed areas of the Konza tallgrass 
prairie. The role of management practices, topographic conditions and 
the uncertainties associated with in situ or airborne surface tempera- 
ture measurements are assessed in this chapter. The second chapter 
deals with evaluation of diurnal and seasonal spectral characteristics 
of the burned and unbumed areas of the prairie. This is accomplished 
based on the analysis of measured spectral reflectance of the grass 
canopies under field conditions, and modelling their spectral behavior 
using a one dimensional radiative transfer model. 



CHAPTER ONE 


Radiative Surface Temperatures of the Burned 
and Unbumed Areas in a Tallgrass Prairie + 


This chapter has been accepted for publication in Remote Sensing of 
Environment. 
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Abstract 


This study was conducted in a natural tallgrass prairie area in the 
Flint Hills of Kansas. Our objective was to evaluate the surface radiative 
temperatures of burned and unbumed treatments of the grassland as a means 
of delineating the areas covered by each treatment. Burning is used to 
remove the senescent vegetation resulting fran the previous year's growth. 
Surface temperatures were obtained in situ and by an airborne scanner. 
These data showed that the unbumed grass canopy was typically warmer than 
the burned one. Measurements of surface energy balance components revealed 
a difference in partitioning of the available energy between the two 
treatments, which resulted in the difference in their measured surface 
temperatures. The magnitude of this difference is dependent on the time of 
measurements and topographic conditions. 
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Introduction 


This study was conducted in preparation for the International Satellite 
Land Surface Climatology Project (ISLSCP) , First Field Experiment (FIFE) , 
which is scheduled for May - November, 1987. The overall objective of FIFE 
is to further evaluate in detail the use of remotely sensed data in 
estimating surface and near surface conditions /properties, based on the 
measured reflected and Quitted radiation, for input to the biospheric 
components of regional and global atmospheric circulation models (Schmugge 
and Seller, 1986). Earlier studies (Carlson, 1981; Price, 1983) have 
evaluated the potential applications of satellite derived surface 
temperatures in estimation of surface energy balance, moisture availability, 
and thermal inertia from surface energy balance models. 

In comparing the accuracy and resolution of ronotely sensed signals 
obtained by airborne sensors with the corresponding ground-measured 
quantities, a complication may arise from their spatial and temporal 
differences. For example, the images that are obtained from airborne or 
space-based systems include primarily the spatial variation in surface/ near 
surface fluxes of reflected or omitted energy, for an instantaneous field of 
view. Characterization of the same areas by comparable ground measurements 
would include both spatial and temporal trends. These differences must be 
resolved for a meaningful interpretation of the remotely sensed signals and 
their effective application for assessing surfaces and near-surface 
conditions/properties . 

The objective of this study was to assess the effects of removal 
(absence) of a layer of senescent vegetation, resulting from growth in the 
previous year(s), on spatial and temporal trends in surface radiative 
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temperatures in a tallgrass prairie. This is one of several requirements 
toward our ultimate objective of developing appropriate methods or models 
that account for the significant sources of variation associated with in 
situ measurements to obtain representative aerial averages for comparison 
with the remotely sensed airborne measurements. 

Materials and Methods 


Site Description : 

The study site was the Konza Prairie Research Natural Area (KPRNA) 
located near Manhattan, Kansas (39° 9'N, 96° 40'W). This area has been 
selected for conducting the ISLSCP-FIFE experiments during May - November 
1987. The predominant soils in this area are silty clay loam with variable 
depth and soil-water-retention characteristics (Bidwell and McBee, 1973) . 

This tallgrass prairie is composed of 3487 hectares of unplowed 
grassland with big bluestem (Andropogon gerardii Vitman) , little bluestem 
( Andropogon scoparius Michx. ) , and Indian grass (Sorghastrum nutans (L. ) 
Nash) as dominant species. Thirty-six other species have been observed in a 
detailed vegetation composition study (The Tallgrass Laboratory, 1984) . 
Since human settlement about 125 years ago, the major treatments on the 
KPRNA have been grazing by wildlife and natural or controlled burning of 
senescent vegetation on some parts of the prairie. 

The area has a temperate, mid-continental climate with warm, moist 
summers and cool, dry winters. Annual mean precipitation is 750 nm, but it 
has a one percent chance of being less than 460 nm or greater than 1400 mm 
(Henderson, 1971) . 

Surface Temperature Measurements 

Surface radiative temperatures of burned and unburned areas of the 
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prairie were measured by handheld infrared thermometers (Everest Intersci- 

ence, Model 110 ) and an airborne scanner (NS001) during summer of 1985. 

The handheld thermometer had a 2° field of view (FOV) and a spectral 

wavelength band of 8-14 yin, and it was held at 45° from nadir at the time of 

measurements. The scanner was on board a C-130 aircraft sponsored by the' 

National Aeronautics and Space Administration, and it had a spectral 

wavelength band of 10.9-12.3 ym a 2.5 milliradian instantaneous field of 

2 

view (100° POV) , which resulted in a pixel size of about 10 m by looking 
straight down (nadir) frcm an altitude of 1500 m. 

Two adjacent areas of the tallgrass prairie were selected for ground 
measurements and extraction of surface temperatures from the airborne 
images. One of the two areas had a layer of senescent vegetation 0.12 m 
thick on the soil surface (unbumed) , which resulted from one or more 
previous years' growth. Green leaf area index of the unbumed area was 
about 1.0. On the adjacent area (burned) , the senescent vegetation was 
removed by means of burning in early spring prior to the resumption of new 
vegetative growth. Green leaf area index of the burned area was about 1.4. 
The two areas were separated by an access road (fire guard) . Vegetation on 
the burned area was more vigorous and less heterogeneous as compared with 
the unbumed area; however, both canopies were about 0.80 m tall. In situ 
surface temperatures were obtained along established transects that covered 
the major sections of the burned and unbumed areas. The surface temper- 
ature data extracted from the aerial images were from pixels located inward 
from the boundaries of the two areas. 


* Tradenames and model number are given for the benefits of reader. NO 
endorsement intended. 
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In situ surface temperature data were recorded by a portable data 
acquisition system (Polycorder Cftmidata , Model 516A) . The handheld infrared 
thermometer was routinely checked before and after each set of measurements 
by measuring the temperature of a reference blackbody target. Air 
temperatures were recorded by a handheld thermometer before, during and 
after each set of surface temperature measurements. 

Components of surface energy balance were measured above the burned and 
unbumed treatments. Net radiation was measured with Swissteco net 
pyrradicmeters . Latent and sensible heat fluxes were measured based on the 
Bowen-ratio surface energy balance approach, using two high-precision 
resistance-thermometer psychrcmeters mounted iirmediately and 1.0 m above the 
grass canopies. Position of the two psychrcmeters were exchanged once every 
6 minutes to cancel the effect of sensor differences in measured fluxes. 
Soil heat flux was measured with heat- flow transducers placed at 0.05 m 
below the soil surface. Soil temperatures were measured, above the soil 
heat-flow transducers, for computing the heat storage in this layer and 
estimating the surface heat flux. Signal from all sensors were recorded by 
Hewlett-Packard data acquisition system composed of an HP-71B microcomputer, 
an HP-3421A analog to digital converter, a HP-9114A disk-drive, and a 
HP-2225B printer. 

Data Analysis 

Digital counts and equivalent calibration temperature data acquired by 
the NS001 scanner on 17-19 June, 1985 were extracted from the computer 
compatible tapes obtained from the Goddard Space Flight Center. Fran the 
calibration data linear regression equations were established between the 
digital counts and the blackbody equivalent temperatures for each scan line. 
Eight sites was selected inside each of the two (burned vs. unbumed) areas, 
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and digital-number values were extracted frcm five contiguous pixels (i.e., a 
total of 40 pixels per treatment area) frcm the image of the scene. The 
blackbody equivalent temperatures computed frcm the digital counts were then 
corrected for emissivity and atmospheric attenuation, using the NESS-80 
atmospheric model (Weinreb and Hill, 1980) . 

The surface energy balance data were processed based on the Bowen-ratio 
approach (Fritschen, 1965) . The procedure suggested by Ohmura (1982) was 
used for filtering the data especially in early morning and evening when 
gradients of sensible and latent heat are very small (i.e. near neutral 
conditions) . 


Results and Discussion 


Error Analysis 

To assess the uncertainties associated with both in situ and aerial 
measurements, we carried out an error analysis. This was deemed necessary 
for quantification of the temperature differences that could be observed 
between the burned and unbumed areas. To compute the uncertainty associ- 
ated with in situ measurements, we considered instrument accuracy and 
precision, surface emissivity, and spatial and temporal variabilities of 
surface temperatures. 

The accuracy and precision of the handheld infrared radiometer provided 
by the manufacturer were + 0.5 °C and + 0.1 °C, respectively. An estimate 
of combined spatial and temporal variabilities was computed frcm the 
available in situ measurements obtained during 17-19 June 1985. A surface 

q 2 2 

temperature variance of 1.5 C was found for a 10 in area, comparable to the 
pixel size of the aerial scanner. The square root of this variance (+1.2°C) 
gives a rough measure of uncertainty for combined spatial and short-term 
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temporal variances associated with the in situ measurements. We assumed the 
emissivity of green vegetation to be 0.98 + 0.01. Typical green plant 
emissivities range fran 0.97 to 0.99 (Oke, 1980). The temperature 
uncertainty associated with +0.01 unit of green vegetation emissivity is 
+0 . 7°C (Phinney and Arp, 1979). The combined root-sum-square (RSS) 
uncertainty associated with in situ measurements by considering all of the 
above factors was found to be +1.5°C. 

Similarly, the uncertainty associated with the factors that affect 
airborne scanner measurements were computed, based on the procedure 
stablished by Cooper et al. (1985) . The instrument precision for the NS001 
scanner was estimated to be +0.2°C, based on the noise equivalent change in 
temperature. Its calibration uncertainty was estimated to be +0.5°C. 
Intra-pixel variability was accounted for by taking the average of 
temperature differences between a central pixel and its eight adjacent 
neighboring pixels. The range of these average values was +0.5°C to +1.0°C. 
We used the conservative value of +1.0°C in our computations . Since the 
airborne scanner was in nadir position, it could view both the vegetation 
and the background (i.e., soil or senescent vegetation). This could result 
in larger variance in the recorded signals. To account for this, we 
assumed a + 0.02 unit change in surface emissivity for the airborne data. 
The uncertainty associated with this emissivity change, based on the earlier 
values, would be +1.4°C. However, if the atmospheric attenuation were 
included in the uncertainty associated with the emissivity, it would be 
reduced to +1.0°C. This is due to a longer pathlength between the sensor 
and the surface, as compared with in situ measurements. The longer 
pathlength includes increased water vapor content in the atmosphere and, 
subsequently, further attenuation of the upwelling thermal energy. This 
results in a reduced effective emissivity of the surface (Price, 1984) . The 
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RSS uncertainty associated with the airborne measurements was estimated to 
be +1.59°C, by taking account of all of the parameters discussed. The data 
used in the error analysis are summarized in Table 1. The aggregate 
RSS uncertainty associated with both in situ and aerial measurements is 
+2.18°C, assuming that these uncertainties are additive. This aggregate 
uncertainty should be used in expressing the aerial scanner temperature data 
in terms of in situ measurements or vice versa. 

The computed RSS value for the NS001 scanner is the theoretical limit 

of uncertainty based on available information about the sensor, and it 

requires some form of validation. The validation was performed using a 

statistical carparision of in situ and scanner data sets. Two coincident 

data sets from 17 and 19 June were used to calculate statistical confidence 

limits derived from a group-comparison t-test. This gives a measure of 

variability of the in situ data as compared with the scanner data. The 

calculated confidence limits for the in situ and scanner data ranged from 

-1.2°C to +2.2°C, which were smaller than or equal in magnitude to the 

aggregate RSS uncertainty. Thus, the variability of the in situ 

measurements would not contribute additional uncertainties to the RSS 

estimate for the scanner. If the confidence limits were larger than the RSS 

value, then RSS estimate would have had to be upwardly revised to reflect 

the additional variability in the in situ measurements. 

In Situ Surface Temperatures 

\ 

In situ, measured, air and surface temperatures for 19 June 1985 are 
presented in Figure 1. During the early hours of the day, surface tempera- 
tures of both burned and unbumed grass canopies were warmer than that of 
the surrounding air. The grass canopy of the burned treatment was warmer 
than that of the unbumed treatment during this period. At midday, both 
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canopies were at equilibrium with the surrounding environment; but, they 
were cooler than the air in the afternoon. The grass canopy of the burned 
treatment was cooler than that of the unbumed canopy during this period. 
This trend was generally observed during the month of June, when soil wa^er 
content was not limiting because of frequent rainfalls (150 mm rain during 
14-29 June) . The observed differences were not statistically significant, 
because of large variance associated with measured surface temperatures, 
especially for the unbumed treatment. 

Similar in situ data collected during the soil-water limiting periods in 
1984 and 1985 showed that grass canopies were warmer than the surrounding 
environments frcm midmorning until midaftemoon (Figures 2 and 3) . The 
surface temperatures of the burned treatment were cooler than those of the 
unbumed treatment during the morning hours , but the trend was reversed 
around midday. This is an opposite trend to the one observed during the 
periods with adequate soil water content (Figure 1) . The standard 
deviations represented with vertical bars in Figures 1-3 depict the extent 
of spatial variability associated with each set of in situ measurements of 
surface temperature along the transects on relatively flat upland areas of 
the tallgrass prairie during a given period. These data suggest a greater 
spatial variation during wet than during dry periods (Figures 1-3) . 

Figures 4 and 5 show the components of surface energy balance measured 
over the burned and unbumed treatments on 23 June 1985. Net radiation flux 

densities for both treatments were comparable (i.e. difference was less than 

_2 

10 V4n ) , partitioning of this available energy, however, was different 
between the two treatments. Sensible heat flux from the unbumed area was 
about 45% greater than the burned area at midday; but, latent heat flux was 
about 18% smaller. Soil heat flux density was also 50% smaller for the 
unbumed as compared with the burned treatment at midday. These differences 
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were primarily due to the layer of senescent vegetation which serves as a 
mulch-layer, over the soil surface in the unbumed treatment. The net 
effect was increased surface radiative temperatures due to reduced evapo- 
transpiration and increased sensible heat flux in the unbumed treatment. 

The standard deviations presented in figures 1-3, however, do not 
include the contributions of other variable features such as topography and 
soil types that exist in the prairie. These features are inherently present 
in data obtained by airborne platforms with large spatial resolutions. 
Figure 6 shows surface temperature differences between the burned and 
unbumed treatments for four major topographical conditions that are present 
in the tallgrass prairie. These include uplands, relatively level land at 
the ridgetops; slopes; toe, a transition zone between slope and lowland; and 
lowland, relatively level land in valleys, but distant frcm stream beds. 
The typical length of the transects established on these four areas was 
300-500 meters. 

The surface temperatures for the unbumed treatment were consistently 
greater than those of the burned treatment, regardless of topography. The 
magnitude of this difference was largest between two treatments for the 
lowland area, followed by toe, slope, and upland, respectively (Figure 6) . 
The circled mean temperatures , however, are not statistically (p=0 . 05) 
significant. The significantly larger surface temperatures difference 
observed between the two treatments in the lowland and toe areas , as 
compared with the upland, can be attributed to less air movement in the 
valley bottom. The soils in the lowland area are typically deep and 
maintain a high level of water content, which results from horizontal and 
vertical seepage of water from the upland and slopes. Thus, lack of soil 
water as a cause of increased surface temperatures in the valley bottom area 
is ruled out; but, shallow soil depth and limited soil water could be the 
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cause of wanner surface temperatures on the slopes as compared with the 
upland area. The air movement on the surface of slopes is considerably 
greater than that in the lowland area. 

Airborne Surface Temperatures 

The surface temperatures obtained from the airborne scanner (NS001) 
measurements also demonstrate the difference between the burned and unbumed 
treatments (Table 2) . These data were obtained during the afternoon hours, 
at about 14:30 CST. The blackbody equivalent surface temperatures of the 
burned treatment were 3.2, 2.6, and 3.0°C cooler than those of the unbumed 
treatment for 17, 18, and 19 June, respectively. These differences 
increased to 3.8, 3.2, and 3.8°C, respectively, after the blackbody 
equivalent temperatures were corrected for emissivity and atmospheric 
attenuation effects, using the NESS-80 atmospheric model (Weinreb and Hill, 
1980) . This is due, in part, to a longer atmospheric pathlength for the 
emitted thermal energy to reach the airborne scanner as compared with the 
handheld thermometer. The increased water vapor in the pathlength results 
in further attenuation of the upwelling thermal energy (Brower et al., 
1976) . These differences are greater than the RSS uncertainty (1.59°C) 
associated with the NS001 scanner. The in situ measured difference between 
the burned and unbumed treatments was approximatly 2.5° C at 14:30 hours 
(Figure 1). This corresponds to a mean difference of 3.8° 2.19°C, based on 
the NS001 scanner data. These data suggest that it is possible to delineate 
the burned from unbumed areas of the tallgrass prairie based on the 
measurements of their emitted thermal radiation by airborne or possibly 
space-based platforms. Burned areas of the grassland can be also 
distinguished from unbumed areas by their spectral reflectance 
characteristics in the visible to middle infrared regions of the spectrum 
(Asrar et al., 1986) . 
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The variance in surface temperatures of the burned treatment was 
consistently lower than that of the unburned treatment, based on both in 
situ and airborne measurements . This is attributed to the layer of 
senescent vegetation that was present on the soil surface in the unbumed 
treatment. This mulch-like layer significantly affects the surface energy 
balance and soil water status of the grass canopy, and also results in a 
more diverse composition of grass versus nongrass species, as compared with 
the burned treatment (Hulbert, 1969; Knapp, 1984) . 

Summary and Conclusions 

This study was conducted to assess the use of surface radiative 
temperatures of burned and unburned areas of a tallgrass prairie as a means 
of delineating the areas covered by each treatment. An analysis of the 
errors associated with the factors that affect both in situ and airborne 
measurements revealed that the uncertainties associated with the data 
obtained by both methods are of the same order of magnitude (+1.5-1. 6°C) . 
This uncertainty is smaller than the existing mean surface temperature 
differences between the two treatments of the tallgrass prairie. The 
magnitude of in situ surface temperature differences between the two 
treatments was dependent on the time of measurements (daily and seasonal) 
and topographic conditions. m situ and airborne surface radiative 
temperatures, however, can be used to distinguish the burned from unbumed 
areas of the tallgrass prairie. These results support our earlier findings 
that the areas covered by each treatment can be distinguished by their 
spectral characteristics in the visible to thermal infrared regions of the 
spectrum. The activities planned for ISLSCP-FIFE will provide more 
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comprehensive data sets for further evaluation of relating in situ point 
measurement to aerial averages. 
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Table 1. Root sum square uncertainties (a) for the handheld radiometer and 
airborne scanner. 


Factor 


Instrument Precision 
Calibration 

Atmospheric Correction 
Spatial Variability 
Location Error 
Bnissivity 


Aggregate Uncertainty 
Root Sum Square 


Handheld Radiometer 
o(°C) a 2 (°C 2 ) 

+0.1 

0.01 

+0.5 

0.25 

N.A. 

N.A. 

+1.2 

1.44 

N.A. 

N.A. 

+0.7 

0.49 


2.19 


+ 1.50 U C 


Airborne Scanner 


a (°C) 

a 2 (°C 2 ) 

+0.2 

0.04 

+0.5 

0.25 

+0.5 

0.25 

N.A. 

N.A. 

+1.0 

1.00 

+1.0 

1.00 


2.54 

+1.59 U C 
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Table 2. Mean black-body equivalent and surface temperatures corrected for 
atmospheric attenuation effect and their associated standard 
deviation for the burned and unbumed treatments of the 
tallgrass prairie. The data were obtained by the NS001 thermal 
scanner on board a C-130 airplane at an altitude of 1500 m in 
1985. 




NO. of 

Surface Temp 

. (°C) 

Date 

Treatment 

Observations 

Black-body 

Corrected 

17 June 

Burned 

40 

24.6+0.75 

27.2+0.87 


Unbumed 

40 

27.9+1.34 

31.0+1.54 

18 June 

Burned 

40 

23.2+0.85 

25.4+1.05 


Unbumed 

40 

25.8+1.35 

28.4+1.67 

19 June 

Burned 

40 

25.7+0.71 

30.9+0.94 


Unbumed 

40 

28.7+1.48 

34.7+2.02 
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List of Figures 


Figure 1. Air and in situ surface radiative temperatures of the burned and 
unturned treatments measured when soil water content was not 
limiting, 19 June 1985. 

Figure 2. Air and in situ surface radiative temperatures of the burned and 
unbumed treatments measured when soil water content was low, 15 
August 1984. 

Figure 3. Same as Figure 2, except for 16 July 1985. 

Figure 4. . Components of the surface energy balance measured over a burned 
flat ridge-top area on 23 June 1985. 

Figure 5. Same as Figure 4, except for an unbumed area. 

Figure 6. Mean in situ surface radiative temperatures of the burned (T B ) 

and unbumed (T y ) treatments and their difference for upland (U) , 
slope (S) , toe (T) , and lowland (L) . Circled temperatures are 
not significantly different, using Fisher's protected least 
significant difference. 
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Abstract 


This study was conducted to evaluate the diurnal and seasonal, 
spectral reflectance characteristics of burned and unbumed areas of a 
tallgrass prairie, based on field measurements and models of radiation 
transport in the grass canopies. Burning of the senescent vegetation, 
resulting from the previous years' growth, is a cannon management 
practice, which results in improved productivity and affects the 
succession of grass species. The burned and unbumed grass canopies 
showed distinctly different, diurnal and seasonal, spectral reflectance 
characteristics in the visible and infrared regions of the spectrum. 
These were attributed to the differences in development of the two plant 
canopies and the azimuthal differences in sensor-sun-canopy positions 
during £ield measurements of spectral reflectance. The radiation 
transfer model properly simulated the diurnal spectral behavior of the 
two canopies. The simulated, seasonal, spectral reflectance values for 
the unbumed grass canopy were greater than the measured ones, because 
of limitations in proper representation of the layer of senescent 
vegetation in the model. 
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Introduction 


Understanding the diurnal and seasonal patterns of spectral 
reflectance in a plant canopy plays a major role in the applications of 
remotely sensed data for inventory and management of natural resources. 
The non-Lambertian properties of plant canopies caused by their archi- 
tecture affect the angular distribution of radiation reflected from 
them. This complicates the development and understanding of 
relationships between attributes of plant canopies and their spectral 
reflectance properties. Our present knowledge of how canopy reflectance 
may vary with different environmental and physiological factors derives 
from a combination of field measurements and modelling efforts. 

A large number of studies have been conducted to investigate the 
effects of surface structure and environmental factors that contribute 
to the observed spectral reflectance characteristics of plant canopies. 
The effects of solar and view angles on plant canopy spectral 
reflectance have been studied for com [1], cotton [2], soybean [3, 4], 
and wheat [5] . These studies showed that nadir and off-nadir reflectance 
measurements for complete or incomplete plant canopies vary greatly with 
the view-zenith and -azimuth angles. The diurnal effects on the 
measured canopy reflectance appeared to be less in the highly reflective 
(near infrared) wavelength bands than in the visible wavelengths. Thus, 
sun-view geometry has a more pronounced effect on the highly absorptive 
red (600-700 nm) wavelength band, especially for incomplete plant 
canopies. Similar results have been reported, based on the modelling of 
radiation transfer in plant canopies [6-11]. 
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Pinter et al. [5] evaluated the effects of plant canopy archi- 
tecture on diurnal and seasonal, spectral characteristics of several 
cultivars of spring wheat. They found that the visible to middle 
infrared reflectance frcm erectophile canopies was usually higher and 
varied more with changing sun azimuth and zenith angles than that from 
planophile canopies. The planophile wheat canopies had a more uniform 
spectral response throughout the day. The effects of solar angle and 
plant-canopy row orientation on spectrally estimated, green leaf area 
index of winter wheat canopies were evaluated by Asrar et al. [12] . 

The objective of this study was to evaluate the diurnal and season- 
al spectral characteristics of burned and unbumed treatments of a 
tallgrass prairie. This is a requirement for proper estimation of 
grass canopy attributes frcm spectral reflectance measurements. 

Materials and Methods 


Site Description : 

Field studies were conducted in 1984 and 1985 in the Konza Prairie 
Research Natural Area (KPRNA) located in the Flint Hills of Kansas (39° 
9'N, 96° 40 'W) . The predominant soil type in this area is silty clay 
loam (Udic Ustall) , typical of the Flint Hills uplands [13] . KPRNA is 
3487 ha of unplowed bluestem prairie, with big bluestem ( Andropogon 
gerardii , Vitman) , little bluestem ( Andropogon scoparius , Michx.), and 
Indian grass ( Sorghastrum nutans (L.), Nash) as the dominant species. 
Thirty-six other species of grasses, forbs, and small shrubs have been 
identified in detailed, vegetation composition studies [14] . 
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The climate of the Flint Hills prairie is temperate , 
mid-continental, with warm, moist sunmers and cool, dry winters. Annual 
mean precipitation is 750 irm, but there is a one percent chance of it 
being less than 460 nm or greater than 1400 itm [15] . 

Two treatments were established in both years, on opposite sides of 
access roads (fireguards) . On one side, the senescent vegetation from 
the previous year(s) was removed by burning in early spring (18-20 
April) , prior to the resumption of new growth. The area on the opposite 
side was left unbumed, resulting in a layer of senescent vegetation 
that covered the soil surface. These two areas will be referred to as 
burned and unbumed, respectively. The two treatments could be accessed 
only by north-south roads, with the burned and unbumed treatments 
located on the west and east sides, respectively, in 1984. The position 
of the two treatments was reversed in 1985. 

Data Acquisition: 

Spectral reflectance and plant canopy parameters were collected in 
1984 and 1985 at the locations described above. Spectral reflectance 

-j- 

measurements were obtained with a Barnes multinodular (model 15-1000) 
radiometer with 15° field of view at several viewing angles (nadir, 15°, 
30° and 45°) , using a truck-mounted platform. This radiometer has three 
discrete wavelength bands in the visible (450-520, 520-600 and 630-690 
nm) , two in the near infrared (760-900 and 1150-1300 nm) , and two in the 
middle infrared (1550-1750 and 2080-2350 nm) regions of the 

+ Trade names are given for reference. No endorsement intended. 
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electromagnetic spectrum. The elevation of the radiometer at the time 
of measurements was 8 m above the soil surface. Canopy reflectance 
measurements were replicated 20 times on each treatment mostly during 
midday (11:00 a.m. to 14:00 p.m. , CST) , and referenced sequentially to a 
painted BaSO^ calibration panel. The spectral measurements were made at 
least once a week (seasonal) and one to four times during selected days 
(diurnal) , when sky conditions were clear. Spectral reflectance and 
transmittance of individual green leaves of grass and non-grass species 
were measured in the discrete wavelength bands that correspond to the 

Barnes radiometer, by a special spectroradicmeter. 

2 

Nine plant samples , each including 0.1 m of ground area, were 
obtained from nine locations on each of the two treatment transects. 
These sites were marked to avoid resampling and measuring of spectral 
reflectance in the future. In the laboratory, the plant samples were 
separated into green grass and non-grass and senescent vegetation. 
Total green leaf area of both grass and non-grass species was 
determined, using an optical area meter (LiCor, model LI-1300) . A cubic 
spline procedure [16] was applied to smooth these measured values of 
green leaf area over time. 


Theory 

The transport of photons through a host medium, such as a plant 
canopy, can be represented by the linearized Boltzmann equation (hereaf- 
ter called the transport equation) . The analytic time-independent 
transport equation is simply a mathematical statement of particle 
conservation and expresses the change in the number of particles in a 
system as the difference between their production (by scattering) and 
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loss (by streaming or by capture) for a point in space at a particular 
energy level. 

We consider a horizontally homogeneous, plane-parallel 
(stratification perpendicular to the Z-coordinate slab-geometry) , 
optically anisotropic, vegetation canopy of finite physical thickness, 
filled densely and randomly with metameric leaf elements . The leaf 

orientation is assumed to be random, i.e., the normals to the upper face 
of the leaves are distributed uniformly along both polar coordinates. 

The slab is considered to be illuminated by collimated beams of 
monoenergetic photons (direct sunlight) and by a diffuse hemispherical 
source (the sky flux) . Further, the slab is treated as being bounded by 
a reflecting soil or senescent vegetation boundary at the bottom (L c : 
the leaf area index of the canopy) . For this slab, we may write the 
following transport equation (a complete derivation is given in [11]) 

y dI(L;£) 

— = I(L;o) - J(L;0) - Q(L;n), (1) 

G dL 

where 

y = cosine of solar zenith angle (9);y e (-1,0) corresponds to 
descending directions and ye (0,1) corresponds to ascending 
directions; 

0. =. unit vector, denoting a direction located in polar coordinates 
by [ y »4» ] , where <l> is the zenith angle; 

L = cumulative leaf area index; accumulated from top downwards; 

I = specific intensity or angular flux. 

The G (£) function determines the area projected on a plane perpen- 
dicular to the direction 0 by a unit leaf area in the canopy [20]. In 
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the case of uniform leaf normal distribution we have 

1 2ir 1 

G(fl) = G = / d+. /dp. | a. -a | =0.5 . (2) 

2ir 0 L 0 L "*■ “ 


The distributed source term is 

iii 2n 1 

J(L;n) = — / d*' / dp' i(L;o' ) p(a' + a) , (3) 

4tt 0 -1 “ 


where w is the single scattering albedo by the leaves. The scattering 
phase function P(£'->a) is 


2 ir 1 

P(n' -a) = 4 / d*. / dp. | cl- -a 1 | I o . a I 

0 L o L -L - 1 ' -L — 1 


and normalizes to 


1 2ir 1 

— / d* / dp p(a' -► a) s l.o, (5) 

4tt 0 -1 ' 


which states that for any incident direction, a scattered photon will 
have unit probability of being deflected into seme outbound path. 

The external first-collision source term Q(L;£) is the sum of four 

independent first-collision source terms, 
to 

Q(L;£) = ttIq PfOQ + a) exp { GL/p Q > + 


a) 2ir 0 

+ "Ih * I dp' p (£‘ *»■ ft) exp { GL/p' } + 

4ir 0 -1 ~ 


to 2ir 1 

+ vl r0 1 d *' f dy ' P ^‘ ■’“£) ex P t - G 0- -L )/p * } + 

4tt 0 0 ^ 
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w Ztt 1 

+ ifl rd f d(j)' f dy' p(n'+ n) exp { -G(L -L)/y* }, (6) 

4ir 0 0 c 

where is a unit vector denoting the direction of incident collimated 
flux. I Q is the intensity of direct sunlight propagating along 8^, and 
it is estimated as 

i 0 = "(sf) / (v } , (7) 

where 6 is the fraction of direct sunlight in the total incident flux 
density P. 1^ is the intensity of diffuse sky flux, assumed here to be 
isotropic and is estimated as 

I d = (l-e)F /it . ( 8 ) 

I is the intensity of isotropically reradiated flux frctn the soil 
surface. I is due to uncollided propagation of direct sunlight 
through gaps in the canopy to the soil surface, 

Vo = -d/*) r s y 0 wI 0 exp { GL c /w 0 } > O) 

where r is the soil hemispherical reflectance. Similarly, I . is due 

S 3Tu 

to uncollided propagation of diffuse skylight to the soil surface, 

0 

Vd * r s 2ir V 7 dy y exp { GL c /y } . (10) 


There is no analytical solution to the transport equation (Eq.l); 
hence, we resort to a numerical solution. The method of discrete 
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ordinates is a very powerful and general tool for the deterministic 
solution of the complicated and refractory transport equation. In this 
method, the angular dependence of the transport equation is approximated 
by discretizing the angular variables p and <J> into a set of NxM 
discrete directions ;4>j> . The set of corresponding weights is 
denoted by [w,;W ^], The spatial variable L is discretized into a set 
of K spatial nodes. The spatial derivative of Eq. 1 is now approximated 
by a finite difference scheme. The cell-centered flux is computed from 
the cell-edged values, using the diamond difference relations. The 
finite difference form of the discrete ordinates transport equation is 
new 


G A 


* I ( L k+l ; -ij) " 1 ^ L k ; — i j ^ * ~ 0,5 * ^ L k+l ; -i.i^ + * “ 


k+Js 


k+l’-ij 


V-ij 


“ j ) ~ QCl-k+yEjj) > k=l,. . . ,K; 1»1,. . . ,N And 

where J(L^ + ^;_^) and Q(L^ + ^/ 2 ; ^j) are discrete analogs of the dis- 
tributed and first-collision source terms. These source terms are 
evaluated by the exact kernel technique, which involves computation of a 
transfer matrix composed of scattering cross sections for every 

transfer (for details see [11]). Equation 11 may be solved 
for I ( L j c+ i/£ij) in terms of 1(1^;^) as 

1 ( L k+1 i j ) = A i 1 ( L k ’—i j ^ " B i { + ^ L k+y-1j) } * y < 0 


and for 1(1^;^) in terms of KL^^^j) as 

1 ( L k *— i j ) = ^ 1/A -j ) 1 ( L k+1 i j ^ + C i { J ( L k+y-ij) + ^ L k+% ; -ij) } ’ 


ji > 0> 


, (ID 


( 12 ) 


(13) 
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where 


A. = { 1 + f i } / { 1 - f i > 


(14) 


B. = { 2f. } / { 1 - f i > ( 15 ) 

C. = { 2f. } / { 1 + f 1 } ( 16 ) 



(17) 


It should be noted here that Eq.12 is used for forward sweep of 
photons into the canopy and Eq. 13 is used for backsweep out of the 
' - canopy. In order to ensure positivity of the ccrnputed flux, the 
negative flux-fixup procedure employed was to simply set the offending 

' - flux to zero. This method was found to be superior to weighting the 

%• „• 

' diamond difference relations [11] . From a knowledge of the boundary 
: Conditions, Eq. 12 and 13 can be used to iteratively solve the transport 
- Equation. A convergence criterion of 10~^ was used for all the 
computation discussed in the following section. Programming details are 
discussed in detail in a previous report [11] . 


Results and Discussion 

Model Validation; 

r. hn Th e discrete ordinates exact kernel techniques for one- and 
r - r two-angle photon transport were coded in FORTRAN- 7 7 as ONEVEG and 
--rde m)V0O r respectively, with the phase function code embedded in them. 
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After testing for photon flux conservation, we did some standard prob- 
lems in atmospheric photon transport to validate the results of our 
codes by comparison with the benchmark results tabulated by Van de Hulst 
[17 J . The results of Van de Hulst are for homogeneous , finite 
atmospheres illuminated by a uniformly syirmetric, conical source with 
vacuum boundary conditions. The spatial variable L in our case corre- 
sponds to the optical depth of Van de Hulst when 0=1, and further, we 
impose vacuum boundary conditions simply by letting r =o. In this case, 
the incident flux density F is equal to unity, and there is no diffuse 
skyflux ( 3 =1) . It should be emphasized here that the results of this 
validation exercise are for the one-angle problem (i.e., the azimuthally 
average of the two-angle problem) . 

Hie phase function used by Van de Hulst [17] is the azimuthal 
average of the Henyey-Greenstein phase function. 


P(y' + y) 


1 2ir {1 " g2) 

2ir 0 d * {l+g 2 -2g | ft. o' |> 3/2 * 


( 18 ) 


where g is the anisotropy parameter. We chose three values for g; 0 
(isotropic scattering), 0.5, and 0.75 (strongly anisotropic for- 
ward-scattering) . Van de Hulst [17] tabulated the solutions of the 
multiple scattering problem as reflection (R) and transmission (T) 
functions. 


R(v>Uq) = 

-I (0 ;y ) / y Q I 0 

; y >0 


(IS) 

T(y .yg) = 

/u 0 I 0 

; y< 0 , 


(20) 


t 


where I (o,y ) is the reflected angular flux at the top and I(L ,y) is 

c 
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diffusively transmitted angular flux at the bottom of the medium. The 

reflection and transmission functions for t = 2 , <0 =1 u = -1.0 and -0.1 

o 

are plotted in Figures 1 through 4. In all cases, the correspondence 
between the results of our discrete ordinates, exact kernel method and 
those tabulated by Van de Hulst is very good; the maximum relative 
difference was less than one- tenth of one percent. 

The reflection functions shown in Figures 1 and 2 display physical- 
ly sound characteristics. The reflected angular flux increases with 
increase in exit zenith angle, and this increase is steeper for oblique 
incident flux zenith angle. This effect arises because, as becomes 

smaller (oblique directions) , the slab thickness increases and, hence, 
more scattering occurs into the reflected direction. However, as the 
slab encountered by the particles becomes thick enough that a signifi- 
cant number of collisions-absorptions occur, the particles disappear 
before reaching the top of the slab (an effect known as limb darkening 
[18]); see Figure 1 for p-*o. For more details on the validation of 
these codes, see our previous report [11]. 

The scattering phase function [P ] of a leaf layer in which 
the leaves are completely at random is plotted in Figure 5. An eight 
order quadrature (DPg) set was used, and only the first 300 points are 
plotted. The phase function is notationally invariant as expected (the 
noise is numerical) . The anisotropy in the scattered flux is evident, 
with relatively more scattered photons introduced in directions (+y ,<{> ) 
close to the incidence direction -* +1] . The phase function is 

symmetric about _£.£ '= o, and, hence, the backscattered fraction is 
equal to the forward-scattered fraction, which is a consequence of the 

to/2, where r L and 

reflectance 


t^ are the leaf hemispherical 


assumption 
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and transmittance. 


Field Measurements; 

The burned and unbumed grass canopies showed distinctly different 
diurnal spectral characteristics in the visible and near to middle 
infrared regions of the spectrum (Figures 5 and 6) . In the visible 
region (630-690 nm) , the measured reflectance of the unbumed treatment 
increased throughout the day for off-nadir viewing angles (i.e., 15°, 
30° and 45°) . This trend was more distinct for the 30° and 45° than the 
15° viewing angles (Figure 6) . The diurnal trend in spectral 
characteristic of the burned treatment was completely opposite. The 
reflectance from the burned canopy decreased consistently throughout the 
day (Figure 7) , except for the rjadir measurement which increased by 1% 
fran early morning to midaftemoon. The decreasing trend was more 
distinct for 30° and 45° than for 15° viewing angles. The diurnal 
change in spectral reflectance of both treatments was least at nadir. 

The opposite trend observed in spectral characteristics of the two 
canopies resulted fron the position of the radiometer with respect to 
the sun and the canopies. The two treatments could be only accessed for 
measurement purposes by a north-south road, with the burned and unbumed 
treatments located on the west and east sides of the road, respectively. 
During the morning hours, the radiometer was positioned in a hemisphere 
opposite to the sun for the burned treatment, recording mostly the 
forward scattered flux of the reflected radiation. The radiometer was 
in the same hemisphere as the sun for the unbumed treatment, and it 
recorded predominantly the backscattered radiation during the morning 
hours. Thus, the azimuthal difference between the sensor-sun-canopy 
resulted in the different trends observed for the visible reflectance 
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characteristics of the two canopies. It should be noted that the 
difference between the canopies was more pronounced in early morning and 
mid-afternoon hours (Figs. 6 and 7) , when contribution from the 
background was minimum, especially for off-nadir measurements. This 
suggests that there were seme basic structural differences between the 
two canopies. Presence of the senescent vegetation on the soil surface 
in unbumed treatment could also have contributed to these differences. 
This is particularly evident in comparing the nadir red reflectance 
(i.e., solid circles in Figs. 6 and 7) of the unbumed canopy which was 
about 3% greater than the burned canopy during the morning hours. This 
difference diminished throughout the day, but the unbumed treatment had 
consistently higher red reflectance for nadir view angle. 

The near infrared diurnal spectral characteristics of the two 
treatments were also different (Figures 5 and 6) . During the morning 
hours, the near infrared reflectance of the unbumed treatment 
decreased, then it increased continuously frem midmoming until midaf- 
temoon. The magnitude of reflectance values for a given viewing angle 
was comparable for the early morning and afternoon measurements. This 
is expected, because of syimnetrical optical properties of the atmosphere 
around solar noon, during clear days. The near infrared reflectance 
from the burned treatment decreased from early morning until midday, 
then it stayed constant until midaftemoon. The measured reflectances 
at nadir showed the smallest change throughout the day for both treat- 
ments . Diurnal spectral characteristics of the two treatments in the 
middle infrared (1550 - 1750 nm) were the same as in the near infrared 
region (Figures 5 and 6) . 

The difference in diurnal characteristics of the two treatments in 
the near to middle infrared regions can be attributed to the azimuthal 
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differences between the sensor-sun-canopy and the layer of senescent 
vegetation that covered the soil surface in the unbumed treatment. 
This layer of senescent vegetation acts as a layer of sponge with large 
interstices, which extinguish the impinging solar radiation reaching the 
bottom of the grass canopy. The separation observed between the near 
and middle infrared reflectance of the unbumed treatment is attributed 
to the opaque character of the senescent vegetation to the exiting 
radiation. 

Measurement vs. Model Results; 

The simulated, diurnal, canopy reflectances of the two treatments 
at an oblique viewing angle (47.9°), for two selected days in 1984 and 
1985, are presented in Figures 7 and 8, respectively. In the visible 
region, the unbumed treatment had consistently higher reflectance than 
the burned treatment for each date, as well as for the two years. The 
difference observed between the burned and unbumed treatments is 
attributed to the layer of senescent vegetation and the difference in 
distribution and magnitude of the green leaf area index of the two 
treatments, especially in 1985. The effect of senescent vegetation in 
the unbumed canopy reflectance is particularly evident in 1984 data, 
since both canopies were at peak greenness. The strong absorption of 
radiation by plant pigments resulted in smaller red reflectance from the 
burned treatment, which had substantially greater green phytanass in 
both years. 

Figure 8 shows the simulated near infrared reflectance fran the two 
grass canopies for 1984 and 1985. In both years, the near infrared 
reflectance from the burned treatment decreased during the morning 
hours, then it plateaued at midday. The reflectance from the unbumed 


45 




treatment showed a consistent decrease throughout the day in both 
years. The unbumed treatment had a higher near infrared reflectance 
than the burned treatment in 1984; however, this trend was reversed in 
1985. This was due to a larger green leaf area (LAI=1.47) for the 
burned than for the unbumed (LAI=0.97) treatment. In this case, higher 
reflectance from the background senescent vegetation was offset by the 
reduced green leaf area of the unbumed treatment in 1985. If the green 
leaf area of the two treatments had been comparable, as in 1984, the 
near infrared reflectance fran the unbumed treatment would have been 
consistently greater than that fran the burned treatment. It should be 
noted that in our modelling effort, we did not include the layered 
structure of the senescent vegetation; it was considered merely as a 
flat surface with uniform spectral characteristics in this region of the 
spectrum. 

The seasonal trends in spectral characteristics of the burned and 
unbumed treatments in the visible and near infrared regions of the 
spectrum are presented in Figures 9 and 10. The connected data points 
are the simulated results from the radiative transfer model and the 
points with standard deviation bars are the measured data at +1 hour 
around solar noon. In general, the red reflectance decreased and near 
infrared reflectance increased with increased green leaf area index of 
the treatments, in both years. This was expected, since development of 
a grass canopy is accompanied by strong absorption of the visible and 
increased reflectance of near infrared radiation. The unbumed treat- 
ment had consistently greater red reflectance than the burned treatment, 
at comparable LAI, both based on the measured and simulated results. 
This was due to the background senescent vegetation, which typically has 
a greater red reflectance than the bare soil [19] . 
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The simulated, near infrared reflectance values for the unbum- 
ed treatment were greater than those of the burned treatment in both 
years /however, the measured values for the unbumed treatment were 
smaller than those of the burned treatment (Figure 10) . This disagree- 
ment between the measured and simulated infrared reflectance data is 
partly due to the method of representation of background senescent 
vegetation in the radiation transfer model. We assumed that the layer 
of senescent vegetation is uniformly flat and Lambertian. This is a 
sinplified representation, since this layer is like a sponge with very 
large interstices, which further extinguish the flux that reaches 
to the bottom of grass canopy. A more realistic representation of this 
layer is currently being explored. 

Sumnary and Conclusions 

The objective of this study was to evaluate the diurnal and season- 
al, spectral characteristics of burned and unbumed areas of a tallgrass 
prairie by direct field measurements and modelling. A one-dimensional, 
photon transport model was developed and tested for standard problems 
with known solutions published in the literature. The model simulated 
accurately the reflection and transmission of photons for all the cases 
considered. 

The measured, canopy spectral reflectance of the burned and 
unbumed grass canopies showed distinctly different diurnal trends in 
the visible and near to middle infrared regions of the spectrum. This 
was attributed to the presence of a layer of senescent vegetation in the 
unbumed grass canopy and the sensor-sun-canopy azimuthal difference 
between the two canopies at the time of measurements. The observed 
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diurnal trends were more pronounced for off-nadir than nadir view 
angles, at all wavelengths for the two canopies. The model simulated 
properly the diurnal spectral characteristics of the burned and unturned 
canopies at an oblique (47.9°) view angle for two different growing 
seasons. 

The measured, seasonal trend in spectral characteristics showed 
consistently greater red reflectance for the unbumed grass canopy than 
for the burned canopy, with comparable green leaf areas. The simulated, 
red reflectance values showed the same trend and compared well with the 
measured values. The measured, seasonal, infrared reflectance of the 
unbumed treatment was smaller than that of the burned treatment, but 
the simulated reflectance values showed an opposite trend. This was 
attributed to a simplified representation of the role of senescent 
vegetation in the simulated results by the model. A more realistic 
representation of the role of this layer is required for proper simu- 
lation of spectral characteristics of the unbumed grass canopy, espe- 
cially in the infrared region. 
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Figure 11. Same as Figure 9, except for the near infrared region of the 


spectrum. 
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